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ABSTRACT 

Have  reflection  by  smooth  media  and  resonance  of  systems  with  radiation 
damping  are  instructive  examples  of  a  failure  of  the  standard  approach  to 
asymptotics*  They  are  also  good  examples  of  a  type  of  exponential  asymptotics 
needed  for  the  sciences*  Successful  modifications  of  conventional*  singular- 
perturbation  theory  have  been  found  for  them  and  show  some  of  the  principles 
promising  wider  usefulness.  They  have  led  to  recent  developments  in  HKB- 
connection  theory,  which  are  also  reported  briefly. 
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SIGNIFICANCE  AND  EXPLANATION 


Asymptotics  plays  a  big  role  in  science  and  technology  because  it  La  a 
form  of  mathematical  analysis  designed  to  produce  highly  illuminating 
results.  The  last  decade  has  brought  a  rapidly  increasing  number  of 
experiences  —  in  studies  of  chemical  reactions  mechanisms,  population 
genetics,  wave  scattering  and  other  matters  of  quite  practical  concern  — 
where  subtle  issues  not  accessible  to  established  asymptotic  methods  turned 
out 'to  be  decisive*  A  few  problems  of  this  kind  have  been  solved,  but  in  a 
way  accessible  only  to  a  small  circle  of  specialists. 

This  article  for  a  book  uses  two  examples  to  explain  to  a  wider, 
mathematically  trained  audience  the  typical  issues  encountered,  the  main  ideas 
which  have  led  to  their  resolution  and  the  changes  in  the  general  approach  to 
asymptotics  which  they  portend. 

One  example  concerns  the  common  process  of  wave  transmission  and 
reflection  in  a  stratified  medium  where,  as  the  wavelength  decreases,  the 
whole  issue  begins  to  slip  unresolved  through  the  net  of  the  established 
mathematical  theory.  The  other,  concerns  unexpected  resonance  effects 
discovered  in  quantum  chemistry,  where  understanding  and  prediction  depend  on 
making  the  old  and  the  new  asymptotics  play  equal  and  complementary  roles. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


have  reflection  and  quasi resonance 


R.  E.  Meyer 


Z.  Wave  Reflection? 

We  are  creatures  of  habit  and  when  a  mathematical  theory  has  enjoyed  an  initial 
success,  it  tends  to  develop  further  under  its  own  steam,  without  much  inquiry  whether  it 
really  addresses  the  questions  most  seriously  at  issue.  It  may  be  of  interest  to  look  here 
at  two  examples  from  mathematical  physics  in  which  it  has  become  apparent  recently  that 
asymptotic  analysis  has  narrowly  missed  the  true  questions  for  a  long  time. 

The  first  example  is  embarrassingly  simple.  When  waves  travel  through  a  medium,  the 
■Min  practical  questions  concern  transmission  and  reflection.  More  often  than  not,  in 
fact,  nothing  else  is  observable.  The  main,  classical  transmission  and  reflection  effects 
arise  from  boundaries,  interfaces,  cracks,  etc.,  where  a  discontinuous  change  in  material 
properties  provides  a  reasonable  model.  Once  such  effects  are  understood,  attention 
wanders  to  the  modulation  of  the  waves  during  travel  through  the  medium,  which  is  not 
normally  a  uniform  one,  but  has  properties  varying  continuously  from  place  to  place.  It  is 
hard  to  suppress  a  feeling  that  those  variations  ere  also  a  plausible  cause  of  partial  wave 
reflection  and  of  correspondingly  incomplete  transmission,  but  little  information  on  it  can 
be  found  in  textbooks. 

To  understand  why,  it  will  help  to  focus  attention  on  the  simplest  circumstances  in 
which  the  salient  points  stand  out  clearly.  Boundaries  and  interfaces  will  be  ignored,  as 
will  be  processes  by  which  one  wave  type,  e.g.  of  compression,  gradually  generates  another, 
e.g.  of  shear.  The  waves  will  be  assumed  linear  and  Fourier-analysed  into  individual  modes 
of  frequency  <■>.  The  material  will  be  assumed  'plane-layared'  so  that  the  phase  velocity 
is  a  continuous  function  c(x)  of  only  one  Cartesian  coordinate  x.  Plane-wave 
propagation  in  the  x-direction  is  then  described  by  as  simple  a  differential  equation  as 
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whore 


d^v/dx^  +  (n(x)/e]^  v(x)  •  0 


<1) 


n(x)  -  Cq/c( x) ,  e  -  Cg/w 

in  tarn  of  a  reference  value  cQ  of  c.  Shortwave  theory  covers,  first  of  all,  the  case 
of  small  wavelength  scale  e ,  but  it  also  tends  to  concern  most  of  what  we  nomally  think 
of  as  waves:  if  e  >  1/3,  say,  the  solutions  will  not  be  recognizably  waves  because  their 
shape  will  change  too  rapidly  fro*  place  to  place  for  an  impression  even  of  the  notion  of 
wave  length.  The  function  n(x)  describes  the  medium  and  has  various  names,  of  which  the 
optical  one,  "index  of  refraction",  will  be  used  here.  Its  definition  makes  it  natural  to 
expect 


n(x)  >  0  (2) 

and  this  will  be  assumed.  For  waves  obliquely  incident  on  the  layered  medium,  the 
following  applies  as  long  as  the  obliqueness  is  not  too  strong  for  ( 2) 1  the  formation  of 
caustics  is  here  excluded  in  the  interest  of  simplicity  (but  it  will  be  central  to  the 
second  example). 

A  classical  definition  of  transmission  and  reflection  also  requires  reference  to 
clear-cut,  unmodulated  wave  states  that  can  be  compared  unambiguously.  This  demands  an 
assumption  that  the  medium  approaches  homogeneity  far  from  the  region  of  notable 
modulation. 


n(x)  ♦  n+  >  0 
♦  n_  >  0 

There  is  then  no  physical  loss  of  generality 
integrable. 


as  x  ♦  “ 

(3) 

as  *♦-*». 

in  assuming  also  that  dn/dx  is  absolutely 


dn/dx  «  L(B)  ,  (4) 

which  is  sufficiently  [Olver  1974]  to  assure  solutions  of  (1)  of  asymptotic,  pure-wave 
character. 


v(x) 


A/15'6 


X 

5  - / 

0 


n(s)ds 
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becomes  large  in  magnitude,  in  turn,  this  justifies  a  radiation  condition# 

v  ^  ^  ^  ^  -#| 

~  t  e1*^®  as  5  ♦  • 

characterising  the  desired  solution  of  (1)  as  an  incident  wave  of  unit  amplitude  plus  a 
reflected  wave  of  amplitude  |r|  on  the  far  left,  and  a  transmitted  wave  of  amplitude 
|t|  on,  but  no  incident  wave  from,  the  far  right.  When  this  condition  is  written  as 

(v  -  ei5/®)ei5/®  ♦  r  as  C  ♦  -- 
-15 /e 


(5) 


ve  ♦  T  as  5  ♦  • 

then  (1)  to  (S)  define  numbers  T  and  r,  the  transmission  and  reflection  coefficients, 
respectively. 

These  two  coaplex  numbers  carry  information  on  both  (real)  amplitude  and  phase,  and 
rather  different  analytical  considerations  attach  to  these  two  aspects.  Questions  relating 
to  phase  will  be  left  aside  here  to  concentrate  attention  on  the  amplitudes  |t|  and 
I r | .  They  are  not  independent,  the  natural  assumption  of  real  index  of  refraction  implicit 
in  (2),  (3)  entails  an  energy-conservation  principle  for  (1)  expressed  by 


\x\*  +  |r*|  -  t  . 


(6) 


The  wave  problem  posed  by  (1)  to  (5)  is  entirely  classical  and  virtually  everything  is 
known  about  its  solution  (Olver  1974] i  it  exists,  is  unique,  and  if  the  limits  (3)  are 
approached  fast  enough,  can  be  described  to  all  orders  by  the  WKB  approximation 


nl/a*  -  ai«/e  l  A  cB  ♦  e'15'®  JBt" 


(7) 


as  e  *  0  for  fixed  S ,  and  by  (2),  the  approximation  is  even  uniform  in  5.  This  ought 
to  furnish  a  reliable  basis  for  the  calculation  of  the  reflection  coefficient,  which  has 
been  carried  out  [Chester  and  Keller  1961]  with  the  following  result. 

WKB-Corollary  1.  If  n(x)  has  k  continuous  derivatives,  except  for  one  finite 
jump  J  of  d^n/dx*  at  xQ,  and  if  d^n/dx**  is  absolutely  integrable  beyond  some  compact 
interval  for  0  <  p  <  k  +  1,  then 

I r |  -  t2n(x0)]“K~1| J|ek  +  oltk)  . 


-3- 


A  brief  proof  is  given  in  the  Appendix*  The  queer  aspect  of  this  result  is  that 
| r |  ia  determined  by  the  jump  of  d'Si/dx'S  regard-leas  of  any  other  properties  of 
n(x),  which  implies  a  further  conclusion  [Schelkunoff  1951]: 

WKB-Corollary  2.  For  a  smooth  index  of  refraction,  with  continuous  and  absolutely 
integrable  derivatives  of  all  orders,  there  is  no  reflection,  |r|  ~  0  to  all  orders  in  e. 

But,  that  is  puzzling  [Mahony  1967]  because  these  theorems  place  no  restriction  on  the 
range  of  n(x),  even  n+  and  n_  need  not  be  close  to  each  other,  and  the  physical 
plausibility  of  partial  reflection  appears  intuitively  more  related  to  the  range  of 
variation  of  the  index  of  refraction  than  to  its  smoothness?  Mathematically,  the  result  is 
equally  paradoxical  because  a  function  in  the  class  C*- 1  can  be  approximated  arbitrarily 

m 

closely  in  any  plausible  norm  by  a  C  -function. 

Generations  have  been  tempted  to  shrug  this  WXB- Paradox  off  as,  perhaps,  merely 
indicating  negligibility  of  reflection  in  smooth  media.  That  will  not  do,  however,  because 
inability  to  calculate  reflection  implies,  by  the  energy  conservation  relation  (6),  that  no 
meaningful  information  on  transmission  is  at  hand  either I 


1(10)  - 0.0005271... 


with  the  third  and  all  further  terms  contributing  less  than  the  last  digit  quoted.  Direct 
computation,  however,  gives 

1(10)  -  -0.0004558...  . 


The  error  of  the  expansion  therefore  exceeds  16%  even  at  n  ■  10,  where  the  expansion  had 
such  excellent  appearance.  Olver  (1964]  points  out  that  this  error  is  closely  accounted 
for  by  the  term 


in  l(n) ,  which  is  technically  negligible  in  comparison  with  all  terms  in  the  expansion, 
but  actually  exceeds  even  l^n4  at  n  «  10. 

In  the  somewhat  larger  context  of  this  Section,  the  WKB-paradox  provides  a  healthy 
cososent  on  a  contemporary  tendency  to  consider  a  problem  solved  when  a  close,  approximate 
solution  of  the  pertinent  differential  equation  and  boundary  conditions  has  been 
obtained.  For  the  simple,  classical  problem  just  described,  we  have  long  known  everything 
about  the  solution,  but  almost  nothing,  about  reflection  and  transisission.  The  solution 
v(x),  however,  cannot  usually  be  observed  inside  the  medium  and  it  has  signally  failed  to 
point  the  way  towards  predicting  what  can  be  observed. 


II.  Central  Scattering 


One  of  the  earliest  and  simplest  problems  of  quantum  mechanics,  which  also  has 

classical  analogues  in  many  sciences,  concerns  the  motion  of  a  particle  in  the  field  of  a 

spherically  symmetrical  potential  U(r).  Its  stationary  states  are  described  [e.g.  Landau 

and  Lifshitz  1974]  by  Schr oedinger ' s  equation, 

*2  2 

t  [E  -  U(r)]T  -  0 

for  the  wave  function  f(x,y,z)e  where  m  is  the  mass  and  B,  the  energy.  It  is 

traditional  to  split  the  angular  momentum  off  by  the  help  of  spherical  harmonics  so 

that  f  ■  r  and  satisfies  a  radial  Schroedinger  equation. 


+  t*  -  0t(r)]* 


0 


(8) 


with  1 centrif ugally  corrected*  potential 

0t(r)  -  u(r )  ♦  fl2t(i  +  1 )/(  2mr  2)  (9) 

where  t  is  the  quantimi  number  of  the  total  angular  momentum. 

A  common  type  of  potential  of  particular  physical  and  chemical  interest  is 
characterised  by  a  central  singularity  of  Ooulomb  type  [Kramers  1926],  so  that 

rtJ(r)  ♦  -0*  <  0  as  r  ♦  0  ,  (10) 

and  by  a  maximum  Da  of  U(r)  at  r  -  rB,  say  (Pig.  1),  whence  U(r)  falls  to  a  finite 

value  as  r  ♦  •,  Which  may  be  chosen  as  D  •  0.  In  physical  parlance,  this  class  of 
potentials  is  defined  by  the  feature  of  a  central  well  surrounded  by  a  potential  barrier 
(Pig.  1).  It  is  well  known  (Landau  and  Lifshitz  1974]  that  bound  states  of  energy  E  <  0 
may  then  exist  in  the  well,  which  are  eigenfunctions  of  (8)  for  eigenvalues  of  E  and 
generate  resonance  in  scattering  processes.  For  positive  energy,  however,  the  effect  of 
tunneling  precludes  bound  states  because  the  leakage  of  probability  through  the  barrier 
Implies  that  any  eigenfunction  would  have  to  decay  in  time.  Indeed,  it  is  not  hard  to 
deduce  rigorously  from  the  quantum  principle  of  conservation  of  total  probability  for 
Schroedinger 'a  equation  that  no  real  eigenvalue  E  >  0  can  exist  for  a  potential  of  this 

type  [Landau  and  Lifshitz  1974]  and  therefore,  no  resonance  can  occur  at  positive  energy. 


In  the  last  decades,  however,  careful  scattering  experiments  have  led  to  the 
observation  of  the  highest  and  sharpest  resonance  spikes  precisely  for  positive  energies I 
The  basic  explanation  of  such  'quasiresonant '  observations  is  not  difficult]  the  leakage 
of  probability  through  the  barrier  implies  radiation  damping  and  therefore,  any  solutions 
of  Schroedlnger's  equation  in  the  tunneling  range  mist  decay  in  time.  In  the  notation  just 
sketched,  they  must  therefore  have  complex  values  of  S,  and  the  characteristic  decay 
time, 

T  -  -ft/ImE  (11) 

is  called  the  life  of  the  solution.  The  tunneling  is  no  one-way  street,  however,  and  as 
outward  tunneling  leads  to  radiation  damping,  so  inward  tunneling  will  produce  a  radiation 
excitation.  The  standard  measure  of  such  excitation  is  called  response  and  is  defined  sb 
follows.  The  generation  of  stationary  states  in  the  tunneling  range  requires  a  supply  of 
radiation  from  infinity  to  compensate  for  the  unavoidable  radiation  damping,  and  the 
'response'  is  the  ratio  of  the  stationary-state  amplitude  in  the  well  to  the  supply 


amplitude  needed  to  maintain  the  stationary  solution.  (Of  course,  these  amplitudes  are 
defined  in  a  mean-square  sense  because  the  local  amplitude  of  ¥  varies  from  point  to 
point.)  It  is  no  great  surprise  to  find,  when  these  notions  are  expressed  more 
quantitatively  [e.g.  Longuet-Higgins  1967]  that  the  response  is  directly  proportional  to 
the  life  T.  In  normal  scattering  experiments ,  accordingly,  solutions  of  short  life  will 
barely  show  up,  but  solutions  of  long  life  may  be  very  strongly  excited.  The  key  problem 
for  the  physicist  and  chemist  is  therefore  the  prediction  of  the  life  (11)  for  those 
solutions  which  have  a  very  long  life. 

Now,  this  tunneling  effect  is  precisely  the  mathematical  problem  for  which  WKB  or 
turning-point  theory  was  first  developed  by  Kramers  [1926],  Zwaan  [1929]  and  banger  [1931] 
and  then  greatly  perfected  by  many  others  [Olver  1974,  1978].  It  achieved  the  formidable 
objective  of  tracing  the  solution  which  is  exponentially  small  in  the  tunnel  reliably 
through  the  shadow  of  the  solution  which  is  there  exponentially  large  in  the  parameter 
A.  The  result  is  an  asymptotic  expansion  of  the  eigencondition  in  powers  of 

k“2  -  A2/<a»Um)  ,  (12) 

whence  an  asymptotic  expansion 

B  ~  l  m0  k”2*  (13) 

s 

of  the  eigenvalues  can  be  deduced,  and  the  problem  appears  solved. 

The  theory  has  indeed  been  successful  for  inelastic  scattering,  both  in  the  quantum- 
mechanical  and  a  classical  context  [Meyer  and  Painter  1979] ,  at  least  as  far  as  the 
determination  of  eQ  in  (13)  is  concerned.  A  closer  look  at  the  fine  print  of  the 
theorems  [Evgrafov  and  Fedoryuk  1966,  Olver  1974]  reveals  technical  difficulties  [Lozano 
and  Meyer  1976],  impeding  the  determination  of  eg  for  s  >  1.  In  any  case,  however,  it 
has  been  shown  [Lozano  and  Meyer  1976,  Meyer  and  Lozano  1983]  that,  for  elastic  scattering. 
every  coefficient  eg  in  (13)  is  real I  It  follows  immediately  from  (13)  and  (11)  that  the 
theory  then  yields  no  Information  at  all  on  Im  B,  on  the  life  T,  on  the  response  and  on 
the  degree  of  quasiresonant  excitation. 


-8- 
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Of  course,  if  the  interpretation  of  this  result  turned  out  to  be,  in  a  manner  similar 
to  the  indications  for  wave  reflection  in  the  preceding  Section,  that  In  E  is 
'transcendentally  small',  Im  E  ~  0  to  all  orders  in  k“2,  then  it  would  be  natural  that 
all  eg  in  (13)  are  real.  It  would  also  mean,  however,  that  the  life  T  and  the  response 
to  excitation  are  transcendentally  Israel  The  eigenvalues  of  imaginary  part  so  negligibly 
small  as  to  fall  through  even  the  fine  meshes  of  turning  point  theory  would  be  precisely 
the  eigenvalues  of  the  greatest  interest. 

If  a  brief  comment  on  the  lessons  of  these  two  examples  be  permitted,  it  appears  that 
the  WXB-expansions  of  the  respective  solutions  v(x)  and  tp(r)  for  them  may  be  the 
correct  answer  to  the  wrong  question? 

The  Author's  experience,  in  fact,  has  been  that  it  is  not  very  rare  that  asymptotic 
expansions  are  of  relatively  little  value  outside  of  mathematics .  That  is  not  real  heresy, 
because  the  basic  concept  of  asymptotics  is  that  of  approximation ,  and  if  the  asymptotic 
property  of  a  first  approximation  be  proven,  then  its  validity  and  value  depend  in  no  way 
on  approximations  of  higher  order. 

From  the  point  of  view  of  mathematical  physics,  the  comment  may  also  be  relevant  that, 
more  often  than  not,  the  solutions  of  the  differential  equations  are  not  themselves  very 
observable.  This  is  canonical  in  quantum  mechanics  and  a  closer  look  at  experiment  and 
field  observation  in  a  number  of  sciences  indicates  that  it  extends  quite  far  into 
classical  physics.  The  main  observables  tend  to  be  quantities  of  the  type  of  scattering 
coefficients  or  resonances,  and  the  two  examples  indicate  that  it  is  not  entirely 
exceptional  to  find  that  their  prediction  requires  approximations  to  both  asymptotic 
quantities  of  algebraic  type  (i.e.,  powers)  and  of  transcendental  type  (e.g., 
exponentials) . 
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Ill .  Have  Reflection 


The  two  examples  just  sketched  were  among  the  earlier  ones  of  an  increasing  number  of 
physical  and  biological  problems  encountered  in  the  last  decade  in  which  asymptotics  of 
exponential  precision  was  found  mandatory.  It  may  therefore  be  of  interest  to  sketch  now 
the  salient  points  of  approaches  that  proved  effective  for  them.  The  analysis  of  wave 
reflection,  in  particular,  has  reached  remarkable  simplicity  and  a  more  general 
significance  of  its  ideas  is  indicated  by  the  surprising  success  with  which  it  has  been 
extended  to  arbitrarily  nonlinear  modulation  [Meyer  1976a]  at  the  instance  of  the  adiabatic 
invariance  of  the  magnetic  moment  in  plasma  physics. 

Of  two  main  steps  by  which  the  reflection  coefficient  |r|  of  Section  I  can  be 
obtained,  the  first  consists  in  no  more  than  the  observation  that  |r|  is  a  number,  which 
must  be  a  functional  of  the  solution  v(x)  of  (1)  to  (5),  and  a  suitable  representation  of 
this  functional  should  be  helpful. 

The  radiation  condition  (5)  indicates  that  the  natural  variable  for  modulation  is  the 
Liouville-Green  variable 


£/e  •  e 


-1 


n(s)ds 


(14) 


which  measures  distance  in  units  of  local  wave  length  and  is  an  analog  of  Hamilton's 
'angle1.  When  the  unknown  v  in  (1)  is  regarded  as  a  function  v(£),  that  equation 

becoates 

d2v/d£2  +  2f ( £ ) dv/d£  +  e-2  v  -  0  , 

1  -2  l15) 
f(?)  -  n  2  dn/dx  , 

and  the  reflection  coefficient  |r|  must  be  a  functional  of  this  modulation  function  f, 
which  is  seen  to  characterize  the  problem  (15),  (5)  completely. 

A  representation  of  that  functional  has  been  obtained  by  many  authors  in  a  variety  of 
ways,  of  which  two  samples  are  quoted  in  [Meyer  1980}.  A  simple  form  of  it  states  that  the 
magnitude,  even  if  not  the  phase,  of  r  is  the  same  as  that  of 


a+  -  /  <[a(C)]2  -  1)e-2i5/ef(5)dC  ,  (16) 

-m 

where  a (i )  is  an  auxiliary  function  defined  by  the  Riccati  equation  of  (1), 

da/dC  -  —  a  +  (a2  -  1)f,  a(-*»)  -  0  .  (17) 

[The  WKB-Corollary  1  of  Section  I  follows  (Appendix)  from  (16),  (17)  by  the  stationary 
phase  rules  for  Fourier  transforms  without  reference  to  the  WKB-representation  of  v  or  to 
details  of  a(C)«]  As  a  pointer  to  the  motivation  for  (16),  it  may  be  noted  that  an 
integral  equation  associated  in  a  simple  and  obvious  way  with  (17)  is 

a(5)e-2i5/C  -  /?  ( [a(s)]  2  -  1)e-2i*/ef(s)ds  .  (IS) 

•jm 

Like  dn/dx,  moreover,  f  (5)  €  L(B),  by  (2),  (3),  (14)  and  (15),  and  what  questions  may 
attach  to  the  integrals  in  (16),  (18)  therefore  tend  to  be  not  questions  of  convergence. 

Since  (17)  indicates  a(€)  to  be  small  —  in  fact,  it  is  readily  proven  to  be  0(e) 

—  a  common  approach  to  the  functional  (16)  has  been  to  iterate  by  the  help  of  (17)  or 
(18),  starting  with  a(€)  -  0,  so  that  a  first  approximation  to  (16)  becomes 

«• 

-  /  e_2i5/ef(5)de  .  (19) 

mm 

That  is  tricky,  however,  because  the  contributions  from  the  oscillatory  integrand  cancel  in 

this  Fourier  integral  with  large  parameter  1/e  to  an  extent  making  the  integral  smaller 

2 

than  the  conjectured  error  0(e  ),  indeed,  smaller  than  any  power  of  e  when  f  (5 )  is 
smooth  and  decays  well  at  " .  The  other  integral  in  ( 16 ) , 

f  a  V«'*f  dt 
2 

turns  out  similarly  to  be  much  smaller  than  o(e  ),  and  in  the  end,  a  correct  execution  of 
this  approach  yields  no  more  than  Corollary  2  of  Section  I,  because  the  functional  (16) 


11- 


possesses  the  favorable  property  of  Fourier  integrals  to  such  an  excessive  degree  that  one 
is  tempted  to  speak  of  cancellation  sickness. 

This  diagnosis  of  the  technical  root  of  the  WKB-paradox  indicates  the  possibility  of 
an  easy  cure  by  the  second  sain  step:  complex  embedding.  To  this  end,  it  is  assumed  now 
that  the  index  of  refraction,  n(x),  is  analytic.  Lest  this  appear  a  drastic  restriction, 
it  may  be  observed  that  n(x)  represents  the  properties  of  the  medium  and  must  be 
specified,  if  not  by  speculation,  then  from  measurements,  which  could  not  support  a 
distinction  between  analytic  and  non-analytic  functions.  A  further  justification  emerges 
from  work  on  a  related  functional  [Meyer  and  Guay  1974,  Stengle  1977]  which  indicates  the 
effective  approach  to  non-analytic,  smooth  functions  n(x)  to  be  their  approximation  by 
analytic  functions. 

When  n(x)  is  analytic,  the  same  follows  from  f(£)  from  (14)  and  (IS),  and  for 
a(E ) ,  from  (17)  or  (18).  A  rational  approach,  in  fact,  is  to  start  from  the  hypothesis 
that  f(C)  is  analytic  on  a  neighborhood  of  the  real  C-axis  of  positive  minimum  width. 

This  demands  an  extension  of  the  radiation  condition  (5)  to  the  analytic  strip  of  f  (E ) »  a 
formulation  is  found  in  [Meyer  1975]  and  permits  shifting  the  path  of  integration  in  (16) 
from  the  real  E-axis  to  a  parallel  path  in  the  lower  half-plane.  On  the  new  path, 

Ira  E  ”  const.  “  -k,  the  offending  factor  exp(-2iE/e)  in  the  integrand  has  very  small 
magnitude  |exp(-2iE/c ) |  -  exp(-2k/e),  and  by  pulling  this  constant  factor  out  of  the 
integral,  the  cancellations  are  made  explicit. 

This  cure  will  be  clearly  improved  by  increasing  k  as  far  as  possible,  i.e.,  by 
shifting  the  path  down  until  it  encounters  the  first  singularity  of  f (E )  (Fig.  2).  For 
simplicity,  only  one  such  singular  point,  ?  »  ?c,  will  be  envisaged  here  (any  finite 
number  of  them  turn  out  [Meyer  1975]  to  make  additive  contributions  to  reflection).  Figure 
2  prompts  a  conjecture  that  a  principle  of  stationary  phase  might  apply  to  the  integral 
(16)  on  this  path,  that  is,  the  contributions  from  the  long,  straight  path  segments  might 
be  of  minor  importance.  This  is  the  first  point  where  the  analysis  calls  for  some  work: 
Simple,  contractive  estimates  on  (18)  [Meyer  1975]  show  the  conjecture  to  hold,  the  major 
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(This  also  explains  why  knowledge  of  the  tnCB-solution  (7)  of  (1)  at  real  x  has  not  been 
relevant  or  helpful  In  the  present  context#) 

The  contribution  of  the  indentation  to  (16)  is  seen  from  (18)  to  be  just  the  junp  of 
a(€)exp(-2i€A )  across  Sc«  so  that  the  remaining  piece  of  the  problem  is  a  local  WKB- 
connection.  Its  solution  is  needed  to  confirm  the  principle  of  stationary  phase  for  (16), 
but  since  the  cancellations  are  already  fully  explicit,  it  is  needed  only  to  a  first 
approximation.  The  solution  [Langer  1931]  has  been  extended  to  a  very  large  class  of 
singularities  in  (Painter  and  Meyer  1982,  Meyer  and  Painter  1983].  The  result  of  the  local 
computation  [Meyer  1975,  1976]  is 
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(20) 


|r|  -  |a+|  »  2  e“^Ccoa(|i*)  +  oja""^6)  , 

whers  ■  »  -la  £  (Fig.  2)  and  the  leae  important  parameter  u  is  related  to  the  branch 
c 

structure  of  f  at  £  .  The  transmission  coefficient  |t|  is  then  given  by  (6). 
c 

The  main  feature  of  reflection  is  now  seen  to  be  the  cancellation  factor  axp(-2a/e ) 
in  (20),  in  which  a/e  aay  be  called  the  'wave  number  characteristic  of  reflection'.  The 
key  paraaeter  a  ■  |ia  $cl  is  ths  halfwidth  of  the  analytic  strip  of  f  (£ )  and,  contrary 
to  intuitive  expectations,  reflection  is  now  seen  not  to  be  closely  related  to  either  the 
range  of  variation  of  the  index  of  refraction  or  to  its  maximal  rate  of  variation.  Though 
clearly  fundamental,  the  width  of  the  analytic  strip  is  a  subtle  property  of  a  function. 

An  interpretation  [Stengle  1977]  that  remains  applicable  well  beyond  the  class  of  analytic 
functions  is  that  a  characterises  the  growth,  as  p  ♦  •,  of  the  L-nora  I  dpn/dxpl  of 
high-order  derivatives  as  function  of  the  order  p  of  differentiation. 

If  the  index  of  refraction  n(x)  be  specified  by  speculation,  a  and  v  in  (20) 
are,  of  course,  readily  read  off  (14)  and  (IS)  (Meyer  1979].  If  the  index  be  obtained  from 
measurement,  however,  the  determination  of  a  to  a  close  approximation  aay  pose  a  problem 
(and  that  of  u ,  aay  thereby  be  aade  moot).  That  this  difficulty  is  peculiar  to  very 
short  wavelengths  is  suggested  by  a  different  approach  to  wave  reflection  [Gray  1982]  which 
accepts  the  restriction  cQf fl  <  1  on  the  modulation  function  in  order  to  solve  (15), 

d2v/d£2  +  e-2  v  -  -2f  (£  )dv/d£ 

for  fixed  e  by  contraction  with  the  simple,  lefthand  resolvent.  In  particular,  a  power 

series  in  ■ fl  usually  provides  an  effectiva  algorithm  for  reflection,  as  long  as  the 

phase  velocity  c(x)  (Section  I)  varies  fairly  slowly.  The  first  tern  in  the  series  for 

2 

(16)  is  then  indeed  (19)  and  the  remainder  is  smaller  by  a  factor  If I  .  Under  normal 
circumstances,  when  the  frequency  is  not  all  that  high  and  the  variation  of  the  index  of 
refraction,  not  exceptional,  reflection  therefore  appears  to  be  more  robust  than  the 
shortwave  result  might  suggest. 


IV.  Quasi resonance 


Tha  central  scattering  problem  of  Section  II  is  technically  harder  and  haa  not  yet 
received  a  treatment  of  comparable  simplicity,  but  a  sketch  of  the  main  notions  and 
principles  by  which  it  was  solved  [Losano  and  Meyer  1976,  Meyer  and  losano  1983]  may  also 
be  of  interest.  The  discussion  of  Section  II  has  served  mainly  to  clarify  that  the 
important,  quasirasonant  states  are  those  of  long  life  (11)  and  that  this  mandates  a  search 
for  eigenvalues  E  of  nonsero,  but  extremely  small,  imaginary  part.  So  small.  Indeed, 
that  it  could  not  be  pinpointed  with  any  conviction  without  rigorous  proof  of  their 
existence. 

Since  the  potential  U(r)  is  real  at  real  radius  r,  it  follows  from  (9)  that  the 
roots  of  E  -  Dj(r)  must  also  be  slightly  complex,  and  since  those  are  the  crucial  turning 
points  of  the  Schroedinger  equation  (8),  it  becomes  clear  that  an  analysis  in  two  complex 
variables,  S  and  r  is  required,  in  combination  with  asymptotics  in  the  real  parameter 
4.  All  experience  to-date  suggests  that  it  may  be  a  principle  of  transcendental-precision 
asymptotics  that  success  depends  on  avoidance  of  prsmature  approximation.  Once  adequate 
conviction  haa  been  attained  that  a  quantity  is  well-defined,  then  it  can  be  given  a  name 
and  the  further  progress  of  the  analysis  need  not  be  impeded  by  the  question  of  how  the 
quantitative  content  of  this  name  might  be  calculated.  Indeed,  it  is  likely  to  become 
clear  only  at  a  quite  advanced  stage  of  the  analysis  which  quantities  really  need  to  be 
computed,  and  to  what  accuracy.  For  quasiresonance,  in  particular,  success  was  first 
achieved  by  conducting  the  analysis  in  the  two  complex  variables  exactly,  if  somewhat 
abstractly,  and  by  postponing  approximation  with  respect  to  11  to  the  very  end.  This  also 
serves  simplicity  by  avoidance  of  entanglement  with  the  details  and  error  estimates  of 
approximation . 

The  first  step  should  clearly  be  to  formulate  the  eigenvalue  problem.  The  governing 
eqution  (8)  can  be  made  non-dimensional  by  measuring  energy  and  potential  in  units  of 

max^lr)  -  UB  -  U(rm) 

(Fig.  1)  and  distance,  in  units  of  rm.  it  then  becomes 

♦  "  +  k2(B  -  Ut)9(r)  -  0,  Ut  -  U(r )  +  1(1  +  1)/(kr)2  , 


(21) 


where  the  large  wave  number  scale  k  is  given  by  (12), 

k  -  ( ^>Um)  1/2rm/«  . 

For  quasiresonance,  attention  may  now  be  restricted  to  angular  momenta  for  which 
2t (t  +  1)/k2  <  max(r3du/dr ) ,  so  that  U^(r)  also  possesses  a  well  (Fig.  1),  to  energies 
in  the  tunneling  range,  0  <  Re  I  <  1,  and  to  wave  functions  satisfying  a  radiation 
condition  that  the  wave  be  purely  outgoing  at  sufficiently  large  |r|. 

Next,  the  potential  0(r)  needs  extension  into  the  complex  plane  of  the  radius  r, 
and  the  reasons  mentioned  in  Section  III  justify  again  a  restriction  to  analytic 
potentials.  More  precisely,  U(r)  is  assumed  analytic  on  an  arbitrarily  narrow 
neighborhood  N  of  (0,*),  beyond  which  it  may  be  left  undefined.  For  a  clear 
formulation  of  the  radiation  condition,  however,  N  is  assumed  'sectorial's  for  all 
sufficiently  large  |r|,  it  is  to  include  an  interval  of  |arg  r|  of  positive  length. 


Figure  3  shows  the  structure  of  the  r- plane  for  smell  In  B  <  0  (as  it  will  turn  out 
to  be,  fortunately,  because  the  analysis  has  no  right  to  assume  its  sign).  There  are  three 
near-real  roots  ra  of  B  -  0 ^ ( r ) «  two  are  plain  from  Fig.  1  and  the  third,  rQ,  lies 
closer  to  the  origin,  where  the  centrifugal  correction  t(i  +  1)/(kr)  over-powers  the 
Coulomb  singularity  (10).  The  figure  also  shows  the  Stokes  lines  of  (21)  on  which 

r  ... 

Re  /  [B  -  Ut(t>]  '  dt  -  0,  s  -  0,1  or  2 
r 

s 

and  of  which  three  issue  from  each  of  the  simple  roots  r8  (except  for  £  »  0,  which  will 
be  ignored  for  a  while,  for  brevity).  The  WXB-theorem  [Evgraf ov  and  Fedor yuk  1966,  Olver 
1974]  associates  with  each  Stokes  line  a  pair  of  fundamental  solutions  u^r),  v^r ) 
of  (21)  which  haws  on  the  character  of  pure  progressive  waves,  undamped  and  un- 

amplified  with  distance  from  r^.  Let  uA  denote  the  wave  outgoing  from  r8  along  L^ 
and  v^,  the  incident  wave.  Both  are  exact  solution*  of  (21)  on  all  of  N,  but  do  not 
possasa  the  pure  wave  character  on  Lj  for  j  *  i. 

The  far-field  Stokes  line  (Fig.  3)  lies  close  to  the  real  axis  and  remains  in 

M,  which  permits  a  precise  formulation  of  the  radiation  condition  that  no  incoming  wave  be 
preeent  at  m i  the  representation 

♦  <r)  -  R.u.lr)  +  B.vjr)  (22) 

of  the  wave  function  as  a  linear  combination  of  u^,  v^  must  satisfy 

B.  ■  0,  M  .  ( 23 ) 

The  final  condition  for  elastic  scattering  is  that  the  wave  function  f  (Section  II)  must 
be  ■  .uare-integrable  and  the  same  follows  for  Y(r).  This  is  effectively  a  regularity 
condition  at  the  singular  point  r  «  0  (Fig.  1),  which  will  emerge  to  be  interpretable  in 
terms  of  the  'reflection  coefficient* 

A0/B„  -  R  (24) 

in  the  central  wave-representation 

*(r)  -  AqUq(t)  +  BqVq ( r )  (25) 


of  the  wave  function 


Since  the  fundamental  pairs  are  exact  solutions,  they  must  be  linearily  related,  and 
since  (22),  (25)  represent  the  same,  exact  solution  4*(r),  it  follows  that  the  amplitude 


coefficients  must  also  be  linearily  related, 

with  a  'scattering*  matrix  S  -  ( <SA j ) )  independent  of  r.  By  (23),  (24),  the  exact 
elgencondltlon  is  therefore 

0  “  ®./®0  “  s22  +  ®21  VB0 

-  S22  ♦  S21  R  . 


(26) 


The  search  for  eigenvalues  is  now  seen  to  involve,  not  the  approximation  of  the  wave 
function,  but  the  'connection'  question  of  how  fundamental  pairs  are  related. 
Schroedinger's  equation  enters  into  observable  matters  only  through  the  three  coefficients 
in  (26),  and  the  only  concern  is  how  those  depend  on  B  and  U(r)  when  k  is  large,  but 
fixed.  The  formulation  chosen  reflects  hindsight  that  this  question  demands  rather 
different  considerations  for  the  singularity  at  r  ■  0,  when  the  angular  momentum  t  is 
bounded  independently  of  k,  and  for  the  scattering  process  away  from  r  •  0,  which  is 
quasiclaasical. 

The  computation  of  the  scattering  matrix  is  precisely  the  objective  of  turning-point 
connection  theory,  which  has  established  several  methods  for  it,  all  leading  to 


V21  ~  1  +  \ 

m 

Y0S22~  •XPt-2X€03{  1  ♦  I  dB(*)k'2} 


(27) 


as  k  ♦  »,  where  YQ  *  0  is  irrelevant  to  (26)  and  is  a  familiar  WKB-distance 

specified  in  (30)  below.  A  definite  algorithm  for  cg  and  dg  has  been  established  only 
under  unrealistic  restrictions  on  the  potential  [Kvgrafov  and  Fedoryuk  1966] ,  but  in  any 
case,  (27)  only  supports  (13)  and  hence,  cannot  yield  any  information  at  all  on  the  life  of 
♦  for  elastic  scattering.  Lozano  and  Meyer  [1976]  therefore  recalculated  S  21  and 
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S22  more  exactly;  since  Olver’s  [1978]  magnificent  'central-connection'  formulation  was 
not  yet  public  at  the  time  of  their  struggle,  they  used  the  'lateral-connection'  approach 
[Bvgrafov  and  Fedoryuk  1966]  in  combination  with  the  principle  of  conservation  of 
probability.  This  makes  E  -  in  (21)  real  when  E  and  r  are  real,  and  permits 

defining  some  of  the  wave  pairs  u^.v^  with  a  complex-con juugate  symmetry  in  the  planes 
of  E  and  r,  which  is  inherited  by  some  of  their  functionals.  E  -  Ut(r)  is  obviously 
analytic  in  E,  moreover,  and  suitable  functionals  inherit  that  also.  By  tracing  this 
analyticity  and  symmetry  painstakingly  through  the  turning-point  analysis,  Lozano  and  Meyer 
[1976]  proved  the  following  result. 

Precision-Scattering  Theorem.  For  potentials  of  the  type  described,  the  scattering 
coefficients  in  (26)  can  be  represented  exactly  in  the  form 


-2k?. 


T0S2i  “  1  expjil^S^/k}  -  (1  +  i){l  +  A(E,k)/k}e  1  , 


(28) 


Tns22  “  *  •xp{iE,<*,k)A}  , 


(29) 


with 


1 

?  <E,k)-/  [F0(r)J1/2  dr  , 

r0 


(30) 


F  -  Ut(r)  -  E 

r2 

?1(B,k)  -  -/  [F^r)]172  dr  ,  (31) 

rl 

(where  the  subscripts  on  F  denote  an  appropriate  determination  of  branches  of  the  root) 
and  with  Yg  *  0,  Z^(E,k)  analytic  in  E  and  real  for  real  E,  and  and  A(E,k) 

bounded  as  k  ♦  •. 

The  crucial,  new  feature  is  here  the  term  in  (28)  with  factor  exp(-2k?1 ),  which  is 
exponentially  small  because  5 1  turns  out  to  be  real  positive  when  E  is  real.  Such  a 
term  would  be  meaningless  in  ( 27 ) ,  but  the  additional  information  on  £  shows  the  first 
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tern  in  YgS^  to  be  of  exactly  unit  magnitude.  The  very  small  term  therefore  describes 
at  real  8  the  difference  between  |YqS21 |  and  unity,  regardless  of  the  much  larger 
uncertainty  about  arg^S^).  It  is  precisely  this  exponentially  small  term  in  (28)  which 
will  emerge  as  the  source  of  all  the  information  on  the  life  T,  and  this  illuminates  how 
the  standard,  technical  meaning  of  larger  and  smaller  in  asymptotics  can  be  misleading. 

So,  it  is  the  principle  of  conservation  of  probability  which  generates  the  symmetry  on 
which  exponential  precision  in  shortwave  scattering  cam  be  founded.  (In  classical 
scattering  (Lozano  and  Meyer  1976],  conservation  of  energy  plays  an  analogous  role.) 

For  short  waves,  the  most  prominent  quantities  in  the  Theorem  are  the  WKB-integrals 
(30),  (31).  An  appropriate  determination  of  branches  of  F1^2  has  been  worked  out  by 
Lozano  and  Meyer  [1976]  and  shows  k£Qexp(-i*/2)  and  k? 1  to  be  real  positive  at  real 
energy  B.  The  former  may  therefore  be  interpreted  as  the  width  of  the  potential  well  of 
Uj(r)  in  (21)  at  the  level  B  in  units  of  local,  radial  wavelength.  Here  not  the 
potential  barrier  just  the  place  where  there  are  no  waves,  kC^  would  be  similarly 
interpretable  as  the  potential-barrier  width  of  in  such  units. 

This  leaves  the  reflection  coefficient  R  in  the  elgencondition  (26)  to  be  analyzed, 
and  Meyer  and  Lozano  [1983]  have  treated  the  case  where  the  angular  momentum  t  is 
'small',  i.e.,  bounded  independently  of  the  wave  number  scale  k,  and  the  case  where  l 
is  so  large  that  IA  +  limit  >0  as  k  ♦  •».  The  latter  case  is  quasiclassical,  the  third 
turning  point,  rQ,  is  a  simple  one,  near  which  the  solution  of  (21)  is  close  to  an  Airy 
function,  the  condition  of  square  integrability  picks  out  the  correct  Airy  function,  and 
the  well-known  solution  (Kramers  1926,  Olver  1974]  of  the  WKB-connection  problem  for  such  a 
simple  turning  point  yields 


R  -  e-i*/2  +  0(k_1)  . 

Unfortunately,  this  result  is  again  inadequate  for  information  on  the  life  T  because 
its  degree  of  accuracy  destroys  the  chance  of  using  the  new  information  of  the  Precision 
Scattering  Theorem  meaningfully  in  the  eigencondition  ( 26 ) .  Nor  would  further  terms  in  the 
asymptotic  expansion  of  R  help  in  that  respect.  But,  Lozano  and  Meyer  [1976]  pointed  out 
that  the  principle  of  conservation  of  probability  for  ( 21 )  permits  normalization  of  the 
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fundamental  pair  Ug.Vg  in  (25)  so  that 

Vr)  “  u0tr) 

at  raal  I,  and  tinea  it  alao  makes  the  wave  function  ♦  real  at  real  B  and  r,  all 
the  unknown  error  tana  must  be  arrangeable  in  complex-conjugate  pairs,  whence  they  deduce 
that  |R|  -  1  exactly  at  raal  energy.  The  analyticity  in  B  then  implies  an  exact 
representation 

R  -  e-t*/2exp[ik_1E0(B,k)]  (32) 

with  Eg  again  analytic  in  B,  bounded  as  k  ♦  •#  and  real  for  real  B. 

The  case  of  small  angular  momentum  is  more  complicated  because  rQ  then  moves  to 
within  0(k~2)  of  the  central,  singular  point  r  -  0  of  (21),  and  that  singularity  now 
over-shadows  the  embryonic  turning-point  structure  at  rQ.  Fortunately,  a  great  deal  is 
known  about  this  Coulomb-singularity  [Kramers  1926,  Olvar  1974]  and  the  connection  results 
for  it  have  been  extended  to  a  large  class  of  other  singularities  by  Painter  and  Meyer 
[1982]  and  Meyer  and  Paintar  [1983].  For  sufficiently  small  |r|,  the  solutions  of  (21) 
are  dose  to  Bessel  functions,  of  which  the  square-integr ability  condition  picks  the 
correct  one,  and  a  sufficiently  careful  comparison  with  (25)  [Kramers  1976,  Meyer  and 
Lozano  1983]  yields 

R  -  e-iW(  ''*+*'  4  (Xk*1'2) 

with 

o<4)  -  4  ♦  ^  -  [4(4  «■  1)] 1/2  .  (33) 

This  is  again  inadequate  for  information  on  the  life,  but  the  same  probability-conservation 
argument  as  for  the  case  of  large  angular  momentum  shows  that  there  must  be  an  exact 
representation 

R  -  e"iW(  1/2  +  2,,exp[ik"1/2E0(B,k)]  (34) 

with  another  function  EQ  of  the  same  properties  as  in  (32). 

This  result  extends  to  4  -  0,  if  the  assumption  [Landau  and  Ufshitc  1974] 
t(0)  *  0  is  added  for  that  case,  and  the  angular-momentum  correction  (33)  to  the  phase 
shift  of  central  reflection  is  then  massive.  For  4*1,  however,  it  is  quite  small  and 
decreases  with  increasing  4,  and  (34)  recovers  (32)  when  t  -  0(k). 
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On*  cannot  help  feeling  that  all  these  technicalities  are  more  complicated  than  they 
ought  to  be  and,  while  hard  analysis  is  surely  unavoidable,  somebody  might  be  able  to 
straighten  it  out  drastically  by  the  help  of  just  the  right  integral  equation  for  (21)?  In 
any  case,  the  results  at  hand  are  sufficient  for  an  exponentially  precise  evaluation  of  the 
eigencondition.  When  (28),  (29)  and  (34)  are  substituted  in  (26),  it  is  natural  to  split 
the  characteristic  form  of  (26), 

A(E,k)  ■  &22  ^  ®21  ^  ' 

into  a  term  collecting  all  the  functions  whence  asymptotic  contribu-tions  of  algebraic  type 
in  k”1  are  to  be  anticipated  and  another  term  that  is  exponentially  small  in  k: 

A(*,k)  -  {AQ<B,k)  +  A1  (B',k)}  ,  (35) 

Aq  -  exp(-2kC0  +  »iO  +  iE^Tc  -  iE0/k1/2]  +  expUE^k)  ,  (36) 

At  -  (i  -  1)(1  ♦  RA)*xp(-2k€1)  .  (37) 

To  establish  now  those  elusive  eigenvalues  I  responsible  for  quasiresonance,  it  is 
convenient  to  begin  with  the  real  roots  Er  of  AQ.  Since  the  appropriate  branch  in  (30) 
makes  €Q(B,k)  «  i |CQ I  for  real  B,  it  follows  straightaway  from  (36)  that  those  roots  are 
given  by 

k|5Q(Br,k)|  +  (Ef  -  t2  +  k1/2EQ)/(2k)  -  (n  ♦  i  +  o)w  (38) 

which  is  just  the  nondimensional  form  of  the  quasiclassical  quantisation  rule  ignoring  the 
radiation  damping  [Kramers  1926,  Keller  1958],  with  Kramers'  [1926]  angular-momentum 
correction  a.  The  new  feature  that  it  is  an  exact  version  of  the  quantization  rule  is  not 
of  much  direct  help,  because  no  practical  algorithm  for  the  evaluation  of  the  E^  has  been 
worked  out.  The  feature  of  immediate  relevance  is  that,  since  U(r)  is  monotone 
increasing  on  (0,rB)  (Fig.  1),  this  quantization  rule  is  known  to  determine  a  unique, 
real  Br(n)  for  large  k  and  given  integer  n  such  that  still  E£(n)  <  Un. 

The  analyticity  in  8  now  permits  application  of  the  principle  of  the  argument 
[ Lozano  and  Meyer  1976]  to  prove  existence  of  a  unique,  simple  root  En  of  A  close  to 
Br(n)  for  all  sufficiently  large  k  and  n.  It  then  follows  immediately  that,  to  a  first 
approximation  as  k  ♦  •, 


(39) 
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•  <* 
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*n  -  *r(n>  ~  -A1(Br,k)/A^(Sr,k) 

-  <1  +  l)tk|5^(Er,k)|)*1«xp[-2k51(Br,k)) 

where 

r, 

ej(«r,k)  -  >1  %it/2  /  |*r  -  Ot(r)r1/2dr  #  0  .  (40) 

ro 

Tha  raal  part  of  (39)  haa  littla  dlract  meaning,  but  tha  Imaginary  part  givaa  the  first 
approximation  to  tha  lifa  (11).  In  tha  original,  dimensional  notation  of  Section  II,  it  is 
T„~  (8m/UB)1/2|5J(*r<n),k)|  exp[ (8mUa)  (Er (n) ,k)/*)  (41) 

which  confirms  tha  conjecture  (Section  II)  that  elastic  scattering  generates  eigenfunctions 
of  a  life  exponentially  large  in  ft.  It  also  shows  that  the  computation  of  such  lives 
requires  no  more  than  evaluation  of  tha  two  definite  integrals  (31),  (40)  of  typical  Mis¬ 
type,  once  the  raal  part,  Er(n),  of  the  eigenvalue  En  has  been  determined  from  the 
quantisation  rule. 


H 


' 


r 


-23- 


h 


The  objective  of  Sections  I  and  IX  was  to  explain  scientific  reasons  for  attention  to 
some  modern  questions  attaching  to  old  and  elementary,  linear  mathematics.  One  problem,  at 
least,  of  fully  non-linear  oscillator  modulation  [Meyer  1976a]  has  greatly  reinforced  those 


reasons.  How  that  initial  answers  to  such  high-precision  questions  have  been  sketched  in 
Sections  III  and  IV,  one  wonders  about  lessons  of  more  general  significance  that  might  be 
drawn  from  them  beyond  those  noted  in  the  preceding  Sections,  namely  relative  unimportance 
of  asymptotic  expansion,  but  importance  of  complex  embedding  and  of  postponement  of 
approximation . 

One  indication  that  has  impressed  the  Author  is  that  the  conventional  comparison 
between  those  asymptotic  contributions  which  are  algebraically  small  and  those,  which  are 
transcendents lly  small,  can  miss  the  point.  Quasiresonance  furnishes  a  particularly  good 
example,  for  the  answer  to  one  of  its  two  key  questions,  viz.  the  eigenfrequency,  depends 
entirely  on  asymptotics  of  algebraic  type,  while  that  to  the  other  key  question,  viz.  the 
resonant  excitation,  depends  entirely  on  asymptotics  of  exponential  type.  There  are 
occasions,  then,  on  which  a  sure  fruitful  view  of  the  distinction  between  'algebraically 
ssull'  and  * transcendentally  small'  may  be  that  this  distinction  is  qualitative  more  than 
quantitative. 

A  second  experience  which  has  Impressed  the  Author  is  that  the  real  observables,  in 
both  examples,  can  be  identified  with  local  properties  of  singular  points  of  the 
differential  equations.  (In  quasiresonance,  most  of  the  points  in  question  are  turning 
points,  but  the  conventional  distinction  between  those  and  singularities  is  all  too 
superficial,  in  any  casei  it  disappears  in  any  intrinsic  formulation,  such  as  (IS).)  It 
would  appear  natural  to  see  a  more  general  significance  in  that  experience,  once  a  complex 
domain  for  the  equations  is  envisaged. 

In  regard  to  wave  modulation  and  scattering,  it  would  also  appear  significant  that  the 
real  concern  of  all  the  hard  analysis,  in  both  examples,  was  not  with  the  approximation  of 
the  solutions  of  the  differential  equations,  but  with  the  connection  of  wave  amplitudes 
across  the  singular  points  of  primary  relevance  to  the  problem.  It  cannot  fail  to  obtrude 


during  the  technical  work,  as  will  surely  have  become  clear  between  the  lines  of  Section 
IV,  that  the  present  form  of  connection  theory  is  laborious,  largely  because  it  involves  so 
much  detail.  The  final  results,  on  the  other  hand,  do  not  really  substantiate  the  need  for 
all  the  detail,  which  has  greatly  discouraged  acquaintance  with  this  branch  of  asymptotics 
and  thereby  made  it  the  preserve  of  a  rather  small  circle  of  specialists.  Does  it  deserve 

the  discredit  or  could  it  be  reformed  to  the  wider  benefit  of  asymptotics? 

Turning-point  theory  is  also  not  very  general,  even  the  great  monograph  [Olver  1974] 
treats  only  the  simplest  types  of  transition  points.  Physics  motivates  such  a  restriction 
in  the  example  of  quasiresonance,  but  not,  in  that  of  wave  reflection.  The  index  of 
refraction  of  a  medium  is  not  ours  to  choose,  but  ours  to  accept  as  we  find  it.  Since  its 

singular  points  dominating  reflection  lie  well  off  the  real  axis  of  distance,  physics 

places  scant  restrictions  on  their  structure.  There  is  no  good  reason  why  they  should 
belong  even  to  the  class  of  'fractional  transition  points'  [Langer  1931,  Olver  1977]. 
Accordingly,  the  mathematical  principle  of  generalisation  might  here  be  helpful  by 
mandating  abandonment  of  detail  and  thereby  promoting  simplicity  and  a  chance  for  guidance 
towards  the  nucleus  of  connection  and  scattering.  Such  an  attempt  has  been  prompted  by  the 
work  sketched  in  Sections  III,  IV,  and  it  may  be  worth  closing  this  article  with  a  brief 
sketch  of  the  results  and  experiences  to  which  it  has  led. 

On  present  evidence,  the  overriding  lesson  seems  to  be  that  wave-amplitude  connection 
may  be  characterized  as  an  asymptotic  expression  of  the  branch  structure  of  the  singular 
point  [Olver  1974,  Meyer  and  Painter  1983]. 

To  carry  this  lesson  from  regular  points  of  differential  equations  [Olver  1974]  beyond 
the  realm  where  detail  is  accessible,  Meyer  and  Painter  [ 1983a]  studied  the  branch 
structure  of  almost  the  whole  class  of  irregular  points  of  linear,  physical  wave-  or 
oscillator-modulation  equations.  In  contrast  to  all  the  earlier  work  on  isolated  singular 
points,  the  new  study  focuses  on  'very  irregular'  points  which  are  branch  points  of 
arbitrary  structure.  The  large  class  of  equations  admitted  is  such  that  each  singular 
point  can  be  linked  by  a  diffeomorphism  to  a  regular  point  of  the  same  differential 
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aquation.  This  lad  then  to  'irregularity  bounds'  on  the  quantitative  degree  of  homotopic 
deformation  of  regular  solution  structure  as  the  diffeoaorphism  is  traced  to  an  irregular 
point. 

An  incidental  discovery  (for  them,  if  not  perhaps  for  every  Reader)  was  that  the 
independent  variable  in  (1)  or  (20)  plays  two  quite  different  roles  in  the  local  solution 
structure  near  the  singular  point.  More  precisely,  this  applies  to  the  natural  variable 
£  in  ( 14 ) ,  which  plays  the  role  of  a  modulation  variable,  while  C/e  plays  the  role  of 
an  oscillation  variable.  Of  course,  this  recalls  immediately  the  notion  of  slow  time  and 
fast  time  in  multiscale  asymptotics.  The  surprise  was  the  discovery  of  it  in  an  analysis 
having  nothing  to  do  with  asymptotical  Meyer  and  Painter ' [1983a]  study  the  'parameter- 
lass'  case  of  the  theory  of  differential  equations,  in  which  e  •  1,  without  loss  of 
generality.  The  two  variables,  moreover,  played  completely  different  roles,  not  in  the 
asymptotic  solution  structure  (which  their  investigation  left  undefined),  but  in  the  local 
structure  at  the  singular  point.  It  would  appear  that  the  multiscale  notion  is  anchored 
much  more  deeply  in  the  singularity-structure  of  a  class  of  differential  equations  than  had 
been  realised  widely. 

The  reason  for  this  foray  into  pure  mathematics  was  the  conjecture  that,  even  in  the 
more  general  context,  connection  is  an  asymptotic  expression  of  local  branch  structure. 
Accordingly,  an  adequate  representation  of  local  structure  should  suffice  for  asymptotic 
connection  of  wave  amplitudes,  and  some  of  the  central  concepts  of  present  turning-point 
theory  might  be  irrelevant  to  that  purpose?  indeed,  the  new  theory  gives  up  both  the  ideas 
of  comparison  equation  and  of  uniform  approximation.  The  reason  is  that  the  class  of 
fractional  transition  points  stands  in  one-one  correspondence  [Langer  1931]  to  the  class  of 
Bessel  functions.  Once  more  general  singular  points  are  admitted,  uniform  approximations 
of  similar  usefulness  cannot  exist.  That  is  a  pity,  for  sure,  but  is  unavoidable  and 
eliminates  temptation  of  detail.  The  comparison  equation  loses  its  usefulness  similarly. 
Instead,  there  is  the  new  idea  of  a  diffeoaorphism  from  regular  to  irregular  points  of  the 
same  differential  equation. 


But,  how  is  asymptotic  connection  to  ba  deduced  from  knowledge  of  no  more  than  local 
structure  at  the  singular  point?  Mayer  and  Painter  [1983]  show  that  the  two-variable 
structure  can  provide  the  key.  They  use  the  'irregularity  bounds'  on  the  extent  of 
departure  of  irregular-point  structure  from  regular-point  structure  to  prove  that  the  two- 
variable  nature  of  the  solutions  assures  distances  from  the  singular  point  at  which  local 
structure  has  not  yet  been  lost,  but  asymptotic  structure  is  already  present.  In  effect,  a 
typical  boundary- layer  concept  has  surfaced  suddenly:  those  bounds  document  'overlap' 
between  the  domains  of  local  and  asymptotic  approximation,  and  it  is  no  great  surprise  that 
the  asymptotic  connection  formulae  then  follow  immediately  from  the  local  branch  structure. 


Appendix 

For  a  brief  proof  of  Cheater  and  Keller's  11961]  WKB-Corollary  1  (Section  1),  it  is 
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again  best  to  ignore  the  MKB-representation  of  the  solution  v(x)  of  ( 1 )  and  to  start  from 
(16),  of  which  a  brief  proof  is  found,  e.g.,  in  [Meyer  1975].  Fron  (14)  and  (15), 

dpf/d5p  -  \  n~P~2  d^’n/dxP*1  +  ...  , 

where  dots  denote  terms  involving  only  derivatives  of  lower  order  than  those  displayed,  so 
that 

[dk”1f/d£lt”1]  -  -1  (n(x„ )  )”'t”'  j 

where  [•]  denotes  the  jump  of  ♦  at  Xg,  [•]«•( xQ  +  <))-♦( xQ-0).  Since 

f(£)  €  L(K) ,  it  follows  from  (17)  or  (16)  that  a(C)  exists  [e.g.,  Ooddington  and 

Levinson  1955],  is  bounded,  in  fact,  is  0(e),  and  has  one  more  continuous  derivative  than 
f(6)  does.  If 

(a2  -  l)f  -  A(5ie)  , 

<sPa/4Zp  -  (a2  -  1)dpf/dCp  ♦  .“  +  2af  dpa/d£p 

-  (a2  -  1)dpf/d5p  ♦  •••  +  2af  a^'lha/e  +  A)/d5p 

by  (17),  so  that 

[dlt"1A/<K,t~1]  -  -  -i  n-k_1  J  (1  +  0(e2)} 

and  dPA/d£p  is  continuous  for  p  <  k  -  2  and  also,  for  p  -  k  -  1  except  at  xQ,  and 
has  absolutely  integrable  skirts  for  p  <  k.  These  properties  support  the  stationary-phase 
evaluation  [Jones  1966]  as  e  ♦  0  of 

e+  -  /  A(£»e)e  ^^d? 

-i m 

to  the  extent  of 

a+  -  (-ie/2)k[dk_1A/dek",J  +  o(ek>  , 

and  | r |  -  |a+|. 
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